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Abstract 

The axial anomaly in the divergence of the singlet axial current in QCD + QED 
^ ■ leads to low-energy theorems for the matrix element of this operator equation over 

{T) . vacuum and two-photon states and for the matrix element over vacuum and two- 

gluon states. The solution of these theorems is related only to the nonperturbative 
phenomena. These matrix elements are calculated in instanton vacuum generated 
N-JL type quark model for arbitrary Nf. It is shown that this model does satisfy the 
low-energy theorems. 



i \ 1. Introduction. 

qj \ The axial anomaly leads to many interesting nonperturbative phenomena in physics. 

^ ■ Among them are 5-violation processes in electroweak (EW) physics, ?7a(1) problem 

in QCD etc. The solution of these problems is intimately related to the topologically 
nontrivial structure of the vacuum in the gauge theories. 

In this paper we apply the axial anomaly low-energy theorems [|J to test the 
chiral quark model |2| which is based on the instanton model of QCD vacuum 0. 
We find that this model does satisfy these theorems. This conclusion provides solid 
background to calculate the different amplitudes of nonperturbative conversion of 
gluons into hadrons and photons. 

2. The low-energy theorems and axial anomaly. 

The axial anomaly in the divergence of the singlet axial current in QCD + QED 
leads to a low-energy theorem for the matrix elements of this operator equation over 
vacuum and two-photon states: 

{0\Nf-£pGG\Zy) = N C ^Y.Q) F{X) F {2) W 

at q 2 = |L|. Here, Nf is the number of the flavors, g is QCD coupling constant with 
2GG = e^ uXo 'G a ^ v G a Xa and G a being the operator of the gluon field strengths , N c is 
the number of the colors, e and Q / are QED coupling constant and the electric charges 
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of the quarks, respectively, 2F^F^ = e^ Xa F^F^J, F$ = e® q iv - e® q^, e^ 2 \ q 1)2 
are polarizations and momenta of photons respectively and q = q\ + q 2 . This relation 
is a consequence of the absence of a massless singlet pseudoscalar boson. Here the 
contribution of the quark masses is neglected. 

So, the problem is reduced to the calculation of the matrix element: 



((% 2 GG|2 7 > = e«4 2 ) J (0\T(g 2 GG J ; m (x 1 )jt m (x 2 )\0) expz^ixj + q 2 x 2 )dx 1 dx 2 . (2) 

As it is evident, gluons can interact with photons only through quark loops. In 
perturbation theory it leads to at least ~ g 4 result for the left hand side of Eq. flU) 
(see e.g. recent discussion of the higher- loop contributions to the axial anomaly ||). 
So, the solution of this theorem is related only to the nonperturbative phenomena 
connected with the structure of QCD vacuum. 

Another nontrivial low-energy theorem concerns the matrix element over vacuum 
and two-gluon states: 

(0\g 2 GG\2gluons) = (3) 

at the limit q 2 = 0. 

These matrix elements are calculated in the instanton vacuum generated N-JL 
type quark model 0, || for arbitrary Nf. 

3. The instanton vacuum of QCD. 

The instanton is the solution of gluodynamics in the Euclidian space 0: 

and in the momentum representation and at small k 2 it is: 



K\k) = OfftJ^^{l + 0{k 2 p 2 )). (5) 



The anti-instanton solution has the same form as in Eqs. (f|) and (||) but with rf m 
instead of ff^ a . Here 0/(1) is the orientation matrix of the instanton I (anti-instanton 
I) in color space, 7/(77) ~ t'Hooft factors 0, p is the size and z is the position of the 
instanton. For large interinstanton distances R » p the sum of the instantons and 
the anti-instantons is also an approximate solution. The calculation of the action for 
gluon fields leads to sum of the actions of free instantons and classical interinstanton 
potential V(R, pi, p 2 , O), where O = 0\0 2 is a matrix of relative orientation. The 
most important part is the instanton-antiinstanton potential Vjj. As it is well-known, 
at large distances it has a form 

32tt 2 Dp 2 p 2 

y n = -r^~ (6) 
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where orientation factor D is 



D 



fi a (l ab n b 



RgRfi 



(7) 



This resembles the dipole-dipole interaction potential and may be attractive. 

The main assumption of the model concerns small distances R ~ p. At small 
distances it is assumed that there is repulsion. This assumption is supported by both 
phenomenological and theoretical consideration ||, || (see recent discussion of this 
"Interacting Instanton Liquid Model" in || and references therein). This leads to a 
stabilization of the size and density of instantons. The distribution of the number of 
pseudoparticles should, for large (N), be given by 



P(N) oc exp 



4 I 8 (N) 



N_ = N + = N/2 



and the effective size distribution of 7's and Ps coincides, and becomes 



d(p) 



const x / 5 exp 



6-4 p 2 



(9) 



When fermions are included (Nf flavors), the coefficient of the beta function 6 is 



6 = — N c - -N f . 
3 3 ; 



(10) 



The quantities (N) and p include effects of the instantons interactions ||. In the 
following calculations, for simplicity, all instanton sizes are considered to be p. 

Both the phenomenological estimates and variational calculations (see recent dis- 
cussion in H ) lead to a mean interinstanton distance of R = (-mxj ~ 1 f m and 
a mean instanton size of p ~ 1/3 fm. The small packing parameter (p/R) 4 = 0.012 
provides a possibility for the independent averaging over positions and orientations 
of instantons. 



4. The chiral quark model. 

The main assumption of this model is the interpolation formula for the quark 
propagator in the single instanton field. It is approximated as the sum of a free 
propagator and an explicit contribution of the zero mode 0, 

(*V(£ /( j))+2m) . - — — ■ 11 

\ y J ' 1— inst 2777, 

Here, V == d — igA, $±(^;^j(j)) is the zero mode wave function of the fermion in 
the background of one 7(7)0. It depends on the collective instanton variables 
- the size p, the position z and the orientation O of the instanton. 
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This interpolating formula should be accurate both at small momenta (p -C 1/p), 
where the zero mode is dominant, and at large momenta (p> 1/p), where the prop- 
agator reduces to the free one. In the background of an -instanton configuration 
and keeping in mind the low density of the instanton media this formula leads to the 
partition function of the model: 



DtfiD^ exp( / (fx ipHdijj) W+ + W-~ , (12) 



where 



/ An 2 p 2 \ Nf fd A z , 
W±= (--/-) J ^det J ± (z), (13) 



v c 



J±U) 



)f9 = J 72^ ex PH(^- 0^)^)^(0^)1(^75)^(0, (14) 

and the contribution of the current quark masses is neglected. The form-factor F is 
related to the zero-mode wave function in momentum space <&±(fc; Cim) and is equal 
to: 

n*) = -4 ftcw*) - hmm - 1 | t - 3 " °_ ^ , (is) 

with £ = \kp. 
The formula: 

(ab) N = J dXexp(Nln^ - N + Xb). (16) 
(here AT >> 1) provides the final expression for the partition function ||: 

Z N = J D^D^exp(-S eff ), (17) 



where 



S eff = J ^idif> + Y+ + Y- 



and 

Y ± = {i) N f\ J d A z det J ± {z) = (-j^) {iMff J d A z det J ± (z). (19) 

The self-consistency condition at the saddle point in Eqs. ([H]), ( |IBD and ( p~9|) leads 
to 

-ftJffjiL (20) 

J (2irYhPF*(k) + B [ ' 



5. Calculations with the low-energy theorems. 

In the quasiclassical(saddle-point) approximation any gluon operator receives its 
main contribution from the instanton background. As an example, for one instanton 
(anti- instanton) 1(1), 
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9 2 G\x) = ro , ^ xol4 = f{x - z), (21) 



[ p 2 + ( x _ z )2] 

and 

S 2 GG(x) = ±/(z-z). (22) 

In the following the operator g 2 GG(x) is considered. Given the low density of the 
instanton medium it is possible to neglect the overlap of the fields of different instan- 
tons. In that case, the matrix element of the gluon operator GG(x) with any other 
quark operator Q is: 



(g 2 GG(x)Q) N = Z^JD^D^exp(J^W) 

x (JV + (W GG+ (x)Q) Wf + -V5" + JV_ (W GG _(x)Q) W^W?- 1 



(23) 



where 



W GG± = ± [~j^-J J d 4 zf(x - z) det J ± (z). (24) 
The application of the formula (|16|) leads to 

(g 2 GG(x)Q) N = J DijD^eMSeff) ((Y GG+ {x) + Y GG _{x)) Q) , (25) 
where 

Y GG± (x) = ± (-) [iM) N f J d 4 zf(x-z) det J±(z), (26) 

and we take N± = N/2 to preserve CP invariance. 

In this particular case, the calculation of the matrix element, Eq. (§), can be 
reduced to the calculation of the the partition function 

Z N [n,a] = Z N l J D^D^ex P (-S eff ), (27) 

with the effective action, S e ff, in the presence of an external electromagnetic field, 
a M , and an external field, n(x), is given as 

- S eff = J ^W^ + Y + + Y_ + Jdx (Y GG+ (x) + Y GG _{x)) k{x), (28) 

where D = d — ieQfd. 
It is clear that 

WWaetoTM-rMW - ^g^_ U.»o (29) 
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Finally, Eq. can be rewritten as: 
-S eff = J^iD^+ 

(f)^" 1 fdz(l+J dxK{x)f{x - z)) det(*MJ+(z))+ (30) 
(if)^' 1 1 dz i 1 ~ I dxK{x)f{x - z)) det(iMJ_(z)). 



6. Bosonization of the partition function Zn[k,o\ 
Another remarkable formula 



exp(A det[L4]) = / dMexp 



-(N f - 1)A ^(detA^)"/" 1 + itr(MA) 



(31) 



is used in the following discussion. It is possible to check this by the saddle point 
approximation of the integral. 
It is convenient to introduce: 

M±(z) = 1 ± / dxn(x)f{x - z)) M 



and rewrite 

-S eff = JtljHDtlj+ 

, Nf-l 



Nt-l 



(f) f ~ ^dzdet(iM + (z)J + (z)) + (^)' J -Jdzd*(iMM)JM)). 
By using the formula, Eq.(RTI), it is easy to show that: 



(32) 



exp((f)' V/ 1 $ dz fet{iM ± {z)J ± {z)) 



J DM± exp 



fdz (-(N f - 1) (f ) 1 (det M ± ) N f-^ +tTr(M ± M ± J i 



(33) 



By using this relation, Eq. (|3"3"|) , the path integral over quarks in the partition function, 
Eq.(|27|), the effective action S e ff, describing mesons in the presence of the external 



fields and k is 



-i 



-S eff [M±,a,K] = fdz (-(JV) - 1) (f) (det M ± ) N f-^ + 

Tr In (iD + iM+MF 2 (1 + («/)) V |(1 + 7s) + iM-MF 2 (1 - (k/))V \{l - 75 

(34) 
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For the processes without mesons the partition function is: 

Z n [k,o\ = 

expTrln (it) + iMF 2 (1 + |(1 + 75) + iMF 2 (1 - («/))*/ * |(1 - 75)) 

^ (35) 
where (k/) = J dxn(x)f(x — z). 

7. The low-energy theorem for the matrix element between vacuum 
and two-photons states. 

The matrix element, Eq. (0), is generated by 

5Z n [k, a] 
5K,(x)5a fi (xi)6a 1/ (x2) K,a ~ 

and is given by the Feynman diagram, Fig.l: 




Fig. 1 

As it is clear from Eq. (53) the factors in the vertices in the diagram are eQf^y^ 
and iMfF^Nj 1 . 

We must calculate A(g 2 ) (in the limg 2 — > 0), which is defined by: 

{2n) 4 5{q -qi- g 2 )A(g 2 ) = / dx exp{-iqx) (0\g 2 GG{x) |2 7 ) = 

e 2 e^4 2) /(0|T(^ 2 GG(x)^ m (xi)^ m (x 2 )|0) expi(-%r + ^xi + q 2 x 2 )dxdx 1 dx 2 . 

(36) 

It is clear from the previous consideration that 

/ 

/■ rf 4 p iMF^sip-^ + tMF^^p + iMF^ip + qz + tMF 2 ) 
X n i (2vr) 4 ((p-g 1 ) 2 + M 2 F 1 4 )(p 2 + M 2 F 3 4 )((p + g 2 ) 2 + M 2 F 2 4 ) 

+ (/i#^,gi# g 2 )]- (37) 
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Here Fi = F((p - qi ), F 2 = F(p + q 2 ), F 3 = F(p) and 

f{(f) = J dxexp(-iqx)f(x) (38) 

is the form-factor of the one-instanton contribution to g 2 GG. At q 2 = 

f(q 2 = 0) = 327T 2 . 
It is easy to show that the trace in Eq.(|37D can be reduced to 

8MV^g 1A g 2CT r(g 2 ), 

where 

rfaW^ ,,. ^,,.« F * „,.„ M> + - (39) 



- (2tt 4 ) ((p - gi ) 2 + M 2 Ff) (p 2 + M 2 F 3 4 ) ((p + q 2 f + M 2 F$, 
It is possible to calculate this integral if we put F = 1. In this approximation 

r( ' > = iSv /„ T^'" 1 + V" ■ (40) 



At a small values of q , we have 



1 f i-*) 



w-towV-mp) (41) 



As a result, the left side of the low-energy theorem, Eq.fll]), is 

w^K^Sa^"'^' («) 

and this coincides with the right side of Eq. ([!]) . 

8. The low-energy theorem for the matrix element between vacuum and 
two-gluons states. 

Here we present, without details, calculations related to Eq. (|3|). This matrix 
element can be written in the form: 



(0\g 2 GG\g(eM, qi ),g(eW,q 2 )} = e g>«eg»* 

(43) 

x Jd 2 d 2 (OlTg^GA^ixJA^x^lO) expi(q lXl + q 2 x 2 )d Xl dx 2 . 

Here A a ^(x) is a total gluon field, q$ are the polarization and the momentum of 
gluons respectively 



As usual, we expand the total field AHx) around the instanton background. The 
main term in Eq . fl4*3|) is the contribution of the instanton background and is order 
~ 0{g~ 2 ). The next term is the contribution of the perturbative fluctuations over 
instanton background and is order ~ 0(g 2 ). It is easy to see from previous consider- 
ations that 0(g~ 2 ) term is given by the formula 

Z N l J D^eMSeff) ((Y GdAA+ (x) + Y GGAA _{x)) Q) , (44) 

where 

Yggaa± = ± (f (iM) N f J d*z f{x - z) 

(45) 

xJdO{-dl)A^{x 1 )(-df)A^(x 2 )detJ ± (z), 

Here the instanton(anti-instanton) is located at the point z with its orientation O. 

Repeating of the bosonization trick leads to the result for the 0(g~ 2 ) contribution 
which is proportional to 

Tr[(id + iMF 2 )~HMF 2 -f S ]. 

It is clear that this Tr and as a consequence 0{g~ 2 ) term are equal to zero. 
The next 0(g 2 ) term is the contribution of the diagrams, Fig. 2. 




Fig. 2 

It is clear that the first diagram is the direct contribution of the operator g 2 GG 
which is equal to 

-g 2 G^G {2 \ 

where 2G^G^ = e^G$ a G ( £ a , Gf u a = efq lu - e®q ifl . 

The factors in the vertices of the second loop-diagram are gX a /2 , ~f^ and iMfF 2, y 5 N 
An account of the contribution from all flavors gives the coefficient Nf. 

A comparison with the previous calculations (Eq.(|36|), (|37|)) leads to the result 
that the contribution of the second loop-diagram is equal in magnitude to the con- 
tribution of the first diagram at q 2 = but of opposite in sign. 

Then, terms of 0(g 2 ) are equal to zero in the limit q 2 — > 0. 

9. Conclusion. 

Thus we conclude that the instanton vacuum generated chiral quark model satis- 
fies the low-energy theorems, Eq.(|]) and @. This provides solid background to cal- 
culate the different amplitude of nonperturbative conversion of gluons into hadrons 
and photons. 
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We are planning to apply this approach to the transitions between the states in 
heavy quarks systems and to the 77 collision processes. 
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